Continuum percolation for Gibbs point 

processes 

Kaspar Stucki*^ 
University of Bern and University of Gottingen 

2 May 2013 

Abstract 

We consider percolation properties of the Boolean model generated 
by a Gibbs point process and balls with deterministic radius. We show 
that for a large class of Gibbs point processes there exists a critical 
activity, such that percolation occurs a.s. above criticality. 

For locally stable Gibbs point processes we show a converse result, 
i.e. they do not percolate a.s. at low activity. 
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1 Introduction 

Let S be a point process in M.'^, d > 2, and fix a i? > 0. Consider the 
random set Zr(E) = IJxes-'^fi(^)' "where Br(x) denotes the open ball with 
radius R around x. Each connected component of is called a cluster. 

We say that S percolates (or R-percolates) if ZrCE) contains with positive 
probability an infinite cluster. In the terminology of [sl this is a Boolean 
percolation model driven by S and deterministic radius distribution R. 

It is well-known that for Poisson processes there exists a critical intensity 
(3c such that a Poisson processes with intensity (3 > f3c percolates a.s. and if 
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13 < 13c there is a.s. no percolation, see e.g. 14|, or [l2| for a more general 
Poisson percolation model. 



For Gibbs point processes the situation is less clear. In [ll| it is shown 
that for some two-dimensional pairwise interacting models percolation occurs 
if the activity parameter is large enough, see also 0] for a similar result 
concerning the Strauss hard core process in two dimensions. 

The aim of this work is to extend these results to any dimension d > 2 and 
to very general Gibbs point processes. We give a percolation condition on 
the conditional intensity of a Gibbs process, which is easily understandable 
and which is satisfied for most Gibbs processes, provided the percolation 
radius R is large enough. Our main result on percolation. Theorem 13.11 
then states, that there exists a critical activity, such that the Gibbs point 
processes percolate a.s. above criticality. The main idea of the proof is the 
contour method from statistical physics. 

Furthermore we state a result on the absence of percolation for locally 
stable Gibbs processes. Here the idea of the proof is that a locally stable 
Gibbs process can be dominated by a Poisson process, and then we use the 
percolation results available for Poisson processes. 

The plan of the paper is as follows. In Section |2] we introduce the Gibbsian 
formalism and give the necessary notations. Our main results are stated in 
Section [3l and in Section Hlthey are applied to pairwise interaction processes. 
Finally, Section |S] contains the proofs of our main results. 



2 Preliminaries 

Our state space is M.'^ for some d > 2 with the Borel a- algebra. Let |A| 
denote the Lebesgue measure of a measurable A C M*^, and let = |IBi(0)| 
be the volume of the unit ball. For two measurable sets A, A' C M.'^ denote 
dist(A, A') = infa:gA,j/eA'||^ ~ vW- If A = {x} we write dist(x. A') instead of 
dist({a;},A'). 

Let (^Tl, J\f) denote the space of locally finite point measures on M*^ equipped 
with the (T-algebra generated by the evaluation maps [^Tl 9 ^ i— )■ ^(A) G Z+] 
for bounded Borel sets A C M'^. A point process is just a D^-valued random 
element. We assume the point processes to be simple, i.e. do not allow multi- 
points. Thus we can use set notation, e.g. x E ^ means that the point x lies 
in the support of the measure ^ and = ClI^'^) denotes the total number of 
points in ^. For a measurable A C M'^ let (OTa^A/a) be the space of locally 
finite point measures restricted to A and its canonical a-algebra, respectively. 
Denote ^|a for the restriction of a point configuration ^ G onto A. 

To define Gibbs processes on R'' we use the Dobrushin-Landford-Ruelle- 
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approach of local specifications, see e.g. j5l,|l0, 13|. Fix a bounded measurable 
A C M*^ and a configuration u G 9Tas called the boundary condition, where 
A'^ = ]R'^\ A. Furthermore fix a measurable function $: ^Tl— T-MUjcxD}, called 
the potential. For a ^ G let 

WA,.(0 = exp [ - (2.1) 
Define a probability measure on by 



ca 



— fl{tj G A}+ V-^ f l{{a:i,...,Xfe}Uw G A} 

\,<^ V fit ^' A A 

ma.o; . . . , Xfc}) (ia;i ■ ■ ■ dxi^ , 



(2.2) 



where 



k=i ' ^ •' ^ 

is called the partition function. Note that the sum in (12. ip and the partition 
function may not exist. However we tacitly assume that the potential $ 
is chosen such that both (12.11) and (12.31) are well-defined for all bounded 



measurable A C M"'; we refer the reader to 13| or for such and related 
questions. A probability measure /i on is called a Gihhs measure, if it 
satisfies the Dobrushin-Landford-Ruelle equation 

/i(A)= / ^A,5|^,(A)^(rfO, (2.4) 

for all A G A/" and for all bounded measurable A C M'^. Let Q{^) denote the 
set of all Gibbs measures corresponding to the potential It may happen 
that Q contains more that one Gibbs measure; such an event is called a phase 
transition. A point process S is a Gihhs point process with potential $ if it 
has a distribution G Q{^). The measures in (12. 2 p are the local specifications 
of /i. They are nothing else than conditional probabilities, i.e. if S ~ /i, then 
we have for all bounded measurable A C M*^, for all boundary conditions 
LO G %A and for a\\ A e M 

P(S G A I S|ac =uj) =^iK,^{A). 



For a; G M"^ and ^ G define the conditional intensity A, see 10 



as 



A(x|0 = exp(-5^<|.({x}ur)). (2.5) 
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Fix a bounded domain A C and a boundary condition u G OTa^^- Then we 
get from (12. ip 

X[x\iUu) = —JT-. , (2.6) 

for all x G A and for all ^ G 9Ta such that x ^ ^ and MA,a;(0 > 0- Let yU G ^($) 
and define TVa = {.^ G 91: ma,c|ac(^|a) = 0}. Obviously fiA,u){N\) = for all 
u G DTa':, and by (12.41) we get fi{N\) = as well. Thus (12. 6p holds for /i-a.e. 
^ G 91 with x^C 

Note that if we are interested in Gibbs point processes on a bounded do- 
main with empty boundary condition (w = 0), Equation (12. 6p coincides, up 
to the null-set A^a, with the definition of the conditional intensity commonly 
used in spatial statistics, see e.g. 0, Def. 6.1]. Roughly speaking, the condi- 
tional intensity is the infinitesimal probability that S has a point at x, given 
that S coincides with the configuration everywhere else. 

For the rest of this paper we assume the potential $ to be constant 
for one-point configurations, i.e. = $({0}) for all x G M*^. Let 

/3 = exp ( - $({0})) and for a; G M'^ and ^ G 91 denote 



A(x|0 = ^^ = exp(- Yl H{^}^0 



The constant (3 is called activity parameter or simply activity. To be able to 
keep track of /3 in our main results, we will mainly use the notation /3A for 
the conditional intensity. 

In order to describe Gibbs processes by the DLR-approach, one can equiv- 



alently characterize them through the conditional intensity, see [10[. There- 
fore, denote Q{/3, A) as the set of Gibbs measures corresponding to the con- 
ditional intensity (3X. 

A Gibbs point process S ~ yU G ^($) is called a pairwise interaction 
process if for every configuration ^ G 9T with |,^| > 3 we have $(0 = 0- By 
denoting Lp(x,y) = exp ( — y})) the conditional intensity simplifies to 

Xix\0=(3ll^{x,y), 

for all X G M'^ and for all ^ G 9t. With a slight abuse of notation, we shall 
use Q{f3,ip) for the set of the corresponding Gibbs measures. 

One of the most important point processes is surely the Poisson process. 
Here, we concentrate only on homogeneous Poisson processes, which can be 
defined as pairwise interaction processes with (p = 1. For Poisson processes 
the parameter (} is called intensity. The more common definition, however. 
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is the following. A point process 11 is called a Poisson process with intensity 
P if for bounded and pairwise disjoint sets Ai, . . . , A„ the random variables 
n(Ai), . . . , n(A„) are independent and Poisson distributed with mean /3|Ai|, 
for i = 1, . . . ,n. It is an easy exercise to show the equivalence of the two 
definitions using fl2.2p . 

For Poisson processes we have the the following percolation result, see 
3, or also 0, Sec. 12.10]. 



Proposition 2.1. Let U be a Poisson process with intensity (3. For all R> 

there exists a critical intensity < /3c < oo such that Zj^{Il) contains a.s. 
only finite clusters if P < Pc o-nd ^^(11) has an infinite cluster a.s. if P > Pc- 



3 Main results 

For our main result on percolation we need the following definition. Let 
S ~ yU G A). We say S satisfies condition (P) with constants r, 5 > if 

A(x I > 5 for all a; G M*^ and for /i-a.e. ^ G 91 with dist(x,^) > r. 

If for all /3 > all Gibbs processes with distribution in Q{P, A) satisfy condi- 
tion (P), then we say that A satisfies condition (P) itself. 

A physical interpretation of condition (P) could be the following. Let 
G 91 and choose a a; G M*^ such that its nearest point in ^ is at a distance at 
least r. Assume that the point process S coincides with ^ everywhere except 
at X. Then the condition (P) states that the energy cost of adding a point 
at X is bounded from above by an universal constant which does not depend 
on the location x nor on the configuration ^. 

The next theorem is our main result on percolation. 

Theorem 3.1. Let A satisfy condition (P) with constants r and 6. Then 
for all R > r there exist a /3+ < oo such that for all P > Pj^ and for all 
Gibbs processes S with distribution in Q{P, A), the set Z^CE) contains a.s. an 
infinite cluster. 

Remark 3.2. In fact, the proof of Theorem 13.11 yields a slightly stronger 
percolation result. Let S be as in Theorem 13. 11 and let H C M"^ be any affine 
subspace of dimension at least two. Then fllHI contains a.s. an infinite 

cluster. 

For our result on non-percolation we need the following stability assump- 
tion. Let S be a Gibbs process with distribution /z and conditional intensity 
A. Then H is called locally stable if there exists a constant c* such that 
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A(x I < c* for all x G and for yU-a.e. ^ E^. If all Gibbs point processes 
corresponding to the conditional intensity A are locally stable, we call A it- 
self locally stable. Most Gibbs point processes considered in spatial statistics 
are locally stable, see P, p. 84] and 0, P- 850]. However the most impor- 
tant example in statistical physics, the Lennard-Jones process is not locally 
stable. 

Theorem 3.3. Let /3A be locally stable. Then for all R > there exists a 
/3- > such that for all (3 < f3- and for all Gibbs processes S with distribution 
in Q{/3,X), the set ZjiCE) contains a.s. only finite clusters. 

Remark 3.4. The proof of Theorem 13 . 3 1 relies on the fact that a locally sta- 
ble Gibbs process can be dominated by a Poisson process, and then Propo- 
sition 12.11 is applied. However for Poisson processes there are more general 
percolation results available, e.g. the balls could be replaced by some random 
geometric objects. Thus Theorem 13.31 may be generalized along the lines of 

The next example combines Theorem 13 . 1 1 and Theorem 13. 3l to get a similar 
result as in Proposition 12.11 

Example 3.5. A Gibbs point process is called an area interaction process, 
see jif, if its conditional intensity is given by 

A(x I = (3^-\^r{x)\Uy^,M,.iy)\ 

for some 7,r > 0. One easily gets the estimates 

1 < A(x I < 7""'"' if < 7 < 1, 

^-a.r" < A(a; I ^) < 1 if 7 > 1. 

Thus A satisfies condition (P) for all > with 5 = min{l, 7""'*^'*} and it is 
locally stable with constant c* = /3 max{l, 7""'*^' }. Then the Theorems 13.11 
and 13.31 yield that for all R> there exists two constants < /3_ < /3+ < 00 
such that all area interaction processes with /3 < do not i?-percolate a.s. 
and all area interaction processes with /3 > /3+ do i?-percolate a.s. It remains 
open, whether /3_ = /J^ as for the Poisson process, see Proposition 12.11 

4 Pairwise interaction processes 

In this section we consider pairwise interaction processes and compare the re- 
sults of pj| with ours. The next proposition shows which pairwise interaction 
processes satisfy condition (P). 
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Proposition 4.1. Letr > 0. Assume that the interaction function (f satisfies 
one of the following conditions. 

(i) For all x,y & with \\x — y\\>r we have (p{x, y) > 1. 
(a) There exist constants 6 > and rmax < oo such that 

(p{x, y) = if < ||x — y|| < r, 

ip{x,y)>6 if r < ||X - y|| < Tmax, 
ip{x,y)>l if Tmax < - 

Then there exists a 5 > such that all Gibbs processes with distribution in 
Q{f3,(p) satisfy condition (P) with constants r and 6. 

Proof. (i) Since X{x \ C,) = Ylye^vi^yV)^ condition (P) is trivially satisfied 
with 6 = 1. 

(ii) Assume 6 < 1, otherwise condition (?) is satisfied. Let fi G Q{f3,(p). 
Fix G M and consider the event 

= G 91: There exist y,y' G ^Ib^^^Jx) with \\y - y'\\ < r}. 

For a bounded A D Br^^^^(x) and any boundary condition u G ^Mr^^^ix^ 
we get fiA,u){Nx) = and by fl2.4p also n{Nx) = 0. Furthermore 
there exists a constant m < oo such that for all ^ G \ we have 
C(®r-max(^)) ^ m does not depend on x (e.g. m can be chosen 

as the maximal number of balls with radius r/2 which can be placed 
in a ball with radius Tmax + r/2). Thus condition (P) is satisfied with 
S = 6"^ for all a; G M°' and for jj-a.e. ^ G 91. 

□ 



In the assumptions on if include our condition (i) of Proposition 14.11 
together with some conditions ensuring the existence of the corresponding 
Gibbs point processes. Furthermore the authors percolation radius R has to 
be greater than \/2r. We could reduce this bound by a factor \/2, but in 
some cases, e.g. for a Strauss hard core process {ip{x,y) = l{\\x — y\\ > r}), 
one would expect a lower bound on the percolation radius of r/2. The main 
difference is however, that the percolation result of [lH is valid only in two 
dimension, whereas our result holds in any dimension d > 2. 

A pairwise interaction process is locally stable in the following two cases. 
Firstly if the interaction function (f{x,y) is bounded by one; such a process 
is called inhibitory. Secondly if it has a hard core radius r, i.e. (p{x,y) = 
whenever ||x — y\\ < r, and (p{x,y) — )■ 1 fast enough as ||x — y\\ — )■ oo. The 
non-percolation result in [ll| treats only the hard core case. However, unless 



for the percolation result, the authors proof is quite different from ours. 
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5 Proofs 



The main idea for the proof of percolation is based on techniques close to 
the contour methqci in lattice models. In particular, it is a modification of 
the arguments in [11] and [l|. Let condition (P) be satisfied with constants 
r, 5 > 0, and let R > r. Choose a m G N such that m > \/d/ {R — r). Divide 
M°' into cubes of length l/m. For this sake define 

Q, = ^^x eR'^ : \\x - 2;||max < 

where we use the maximum norm ||x||niax = niaxi=i^...^rf|xj|. Then {Qz, z G 
— Z'^} covers the whole space R''. Two cubes Qz and Qz' are called neighbours 
if II -s — 2;' II max = l/m. A set S C —U^ is called connected if for each pair 
{z, 2;'} C S there exists a sequence of neighbouring cubes Qzi, Qz-^-, ■ ■ ■ , Qz^ 
with z = zi and z' = Zn- The cardinality of 5" is denoted by 15*1. 

Lemma 5.1. There exists a constant c > such that for all S C — Z'^ with 
\S\ < 00, there exists a S' G S with \S'\ > c\S\ and for each pair z,z' G S' 
we have dist{Qz,Qz') > r. 

Proof. We can assume 15*1 > 1 and then obviously there exists a 5" with 
cardinality one. To identify the constant c, consider the following inductive 
procedure. Choose an arbitrary point zi & S and draw a ball around zi of 
radius r + \/d/m. Exclude all points of S which are contained in this ball 
and denote by Si the remaining points of S. Note that dist^Qz^, Qz) > ^ for 
all z E Si. Continue this procedure and define S' = {zi, . . . , Zn}, the set of 
the chosen points. At each step we exclude at most ad{r + 3\/d/{2m)Ym^ 
points. Thus 

1^1 



n > max 



d 



where [xj denotes the greatest integer less than or equal to x. Choosing 
c = l/{2ad{r + 3\/d/ {2m)Ym'^) yields the claim. □ 



The following Lemma is the key ingredient for the proof of Theorem 13. 1[ 

Lemma 5.2. Assume S fi E Q{P, X) satisfies condition (P) with constants 
r and 6. Let S C — Z"' with \S\ < 00. Then if P > we have, 



p(dist(H,U,5^?.)>0< (^) 



-c\S\ 

(5.1) 



where c is the same constant as in Lemmal5A[ 
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Proof. By Lemma l5.ll choose a 5" C S* with \S'\ > c\S\ and such that for 
each pair z,z' G S' we have dist{Qz, Qz') > r. Set n = \S'\ and denote 
S' = {zi, . . . , Zn}- Consider the events 

As,r = {^em : dist{^,[j^^sQ^)>^} and 

Bs,r = {^U {xi,...,Xn} e Dft: ^ G As,r, Xi G Qz,, for i = l,...,n}. 

Choose a bounded A C M*^, such that the distance between S and the bound- 
ary of A is at least r + y/d/m and fix a boundary condition w G OTac Then 
for /i-a.e. ^ U {xi, . . . , x„} G -85.,. we have by fl2.6p and condition (P) that 

MA,a;(^|A U {Xi, . . . ,X„}) = Ma.o; (^|a) A (xi | U w) X 

A(X2 I ^|a U {Xi} Uw) ■ ■ ■ A(x„ I ^|a U {Xi, . . . ,X„_i} U Uj) 

> WSruA^A). (5.2) 
The partition function can then be bounded from befow as 

CA.CO > 77 / ■ ■ ■ / ^{{Xl, . . . , Xfc} G Bs,r} X 

uk,ui{{xi, . . . ,Xfe}) dxi ■■■dxk 



X 



. . ,Xfc} G As^r]uK^^{{xi,...,Xk}) dxi---dxk 

ua,uj{{xi, . . . ,Xj}) dxi ■ --dxj, 

where k\/{k — n)\ is the number of possibihties to choose a set of A; — n 
elements and n sets of one elements out of a set of k elements, and the last 
inequahty follows by (15. 2p . by integrating over the cubes {Qz,,^ = 1, . . . , n} 
and by the change of variable j = k — n. Thus, by (12. 2p 

f.AAMr)<{^y^. (5.3) 

Since (15. 3 p does not depend on u, Equation (12. 4 p yields 

P(d,st(= U,.Q.) > r) = P(E e M,.) < < (^)'"". 

where the last step follows hj n > c\S\. □ 
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Remainder of the proof of Theorem \3.1[ We call a set L C — Z'^ a loop if 
\L\ < oo, L is connected and each cube {Qz, z E L} has exactly two neigh- 
bours. We say L C ^Z'^ is a loop around the origin if the origin is contained 
in the convex hull of L. Denote by Cq the set of all loops around the origin. 

Consider the loop L = {zi, . . . , Zk} G Co such that Qz^ and Qz,, are neigh- 
bours and for i = 1, . . . , A; — 1 the cubes Qz^ and Qz^+i are neighbours. Since 
the origin is contained in the convex hull of L, the point zi lies necessar- 
ily in the big cube [—k/m,k/m]'^, which gives {2k + lY possibilities. For 
i = 1, . . . , k — 1, the cube Qzi+i is a neighbour of Q^., thus there are at most 
2"^ possibilities to place the cube Qzi+i- We conclude, that there are at most 
{2k + iY2'^^''~^^ loops in Cq with length k. Consider the event 

Ak = {C, E D^I: There exists a L E Cq with \L\ = k and dist(^, [Jz^lQ^] ^ ''^}- 
Then for /3 > Lemma [5.21 yields 

P(HGAfc)=p( U {dist(S,U,,W.) >r}) 

LeCo,\L\=k 

< Y: F(dist(H,U,,g.) >r) < (2^ + l)^2'^(-^)(^)~". 

LeCo, \L\=k 

Set /3+ = {2^/^mY/6 and choose a (3 > f3+. Then, 

°° °° / Rr)\ -ck 

J2 m e A,) < ^(2fc + l)'^2'^('=-i) (^) < oo. (5.4) 

k=l k=l 

The first Borel-Cantelli Lemma together with fl5.4p then yield that there are 
a.s. only finitely many loops L E Cq with dist(S, IJ^gL Qz) > "r- 

The following trick allows us to reduce the problem to the two-dimensional 
case. Consider the set Z^^\~) = Zr{E) n R^, where we identify as the 
plane spanned by the first two canonical basis vectors in W^. A loop is said 
to be a M.'^-loop if the centres of its cubes lie in M^. 

The parameter m is chosen such that if two points x and x' lie in the 
same connected component of \ Z^^^(S), then there exists a sequence of 
neighbouring cubes with centres S = {zi, . . . , Zn} C ^Z'^flM^, with x G Qz^, 
x' G Qz„ and dist(H, IJ^^^ Qz) > r. To see that this is possible consider a 

path in \ which connects x with x' and choose all cubes which 

intersect with this path. Obviously the cubes are connected and the length 
of the diagonal of a cube is less than R — r. Thus, by the reverse triangle 
inequality, the selected cubes and S are separated by a distance at least r. 
Figure [1] shows a graphical illustration of this procedure. 
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Figure 1: Here we choose d = 2, r = 0.2, R = 0.3 and m = 15. The grey area 
corresponds to and the chain of shaded boxes shows one possibihty to 

separate the two components of Note that the boxes do not intersect 

with the circles of radius r. 

Assume that Zr(E) does not percolate, i.e. there exists a.s. only finite 

(2) 

clusters. Obviously the clusters in Z)j (S) are also finite, and each cluster 
can be surrounded by a R^-loop L E Cq with dist (S, IJzgL Q^) — ^- large 
/3 this eventually leads to a contradiction. 

Note that the choice of a basis in is by no means important. Thus, 



Proof of Theorem \3.3[ Let S be a locally stable Gibbs process with condi- 
tional intensity /3A < c*. Denote A„ = [— n, n]'^, for n G N. Note that 
conditioned on S|a^ = uj, the law of S|a,^ is uniquely characterized by the 
conditional intensity /3A(- | ■ Uw). 

It is a known fact that every locally stable Gibbs process on a bounded 
domain can be obtained as a dependent random thinning of a Poisson process; 
see 0, Remark 3.4]. In particular, there exists a Poisson process n„ with 




□ 
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intensity c* and a Gibbs process H„ ~ G{f3X{- \ ■ U u)) on A„, such that 
S„ C n„ a.s. Let S„ be the point process obtained replacing S|a„ with H„, 
i.e. H„ = S|ac UH„. Since H|a„ and S„ have the same conditional distribution 
/^An.w, Equation (12.41) yields that S and S„ have the same distribution for all 
n e N. A standard result ([3|, Theorem 11.1. VII]) then yields that as n — )■ cxd 
the sequence S„ has a weak limit S' with the same distribution as S. By the 
same result n„ — > 11 as n — t- oo, where 11 is a Poisson process on M°' with 
intensity c*. Furthermore, since S„ C n„ a.s. for all n G N, we conclude 
S' C n a.s. 

Obviously Z^iiE') C 2^(11). Thus if Zuill) does not percolate, neither 
does Zji(E'), and Proposition 12. II finishes the proof. □ 
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